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Abstract 

A search for heavy Majorana neutrinos in events containing a pair of high-px leptons of the 
same charge and high-px jets is presented. The search uses 20.3 fb“^ of pp collision data 
collected with the ATLAS detector at the CERN Large Hadron Collider with a centre-of- 
mass energy of -\A = 8 TeV. The data are found to be consistent with the background-only 
hypothesis based on the Standard Model expectation. In the context of a Type-I seesaw 
mechanism, limits are set on the production cross-section times branching ratio for produc¬ 
tion of heavy Majorana neutrinos in the mass range between 100 and 500 GeV. The limits 
are subsequently interpreted as limits on the mixing between the heavy Majorana neutrinos 
and the Standard Model neutrinos. In the context of a left-right symmetric model, limits on 
the production cross-section times branching ratio are set with respect to the masses of heavy 
Majorana neutrinos and heavy gauge bosons Wr and Z'. 
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1 Introduction 

The discovery of mixing between generations of neutrinos [1] has established that at least two of the 
neutrinos have small non-zero masses. A unique feature of neutrinos compared to other fermions in the 
Standard Model (SM) is that neutrinos could be their own anti-particles, so-called Majorana fermions. If 
this is realised in nature, then the unusually low mass scale of the light neutrinos could be generated by 
a seesaw mechanism [2-7] which would imply the existence of yet unobserved heavy Majorana neutrino 
states. The nature of Majorana neutrinos would explicitly allow for lepton number violation. 

In this paper, a search is presented for heavy Majorana neutrinos using the ATLAS detector at the Large 
Hadron Collider (LHC). The data sample was collected in 2012 during TeV pp collisions and 

corresponds to an integrated luminosity of 20.3 fb“'. Heavy Majorana neutrinos with masses above 
50 GeV are considered. In this regime, the production and subsequent decay of heavy Majorana neutrinos 
could lead to a final sfafe confaining exacfly fwo charged leptons, where fhe lepfons may have fhe same 
or opposife charge in equal fracfions of fhe heavy neufrino decays. Only lepfon pairs of fhe same charge 
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(same-sign) are considered as there is a smaller expected SM background for pairs of same-sign leptons 
than for pairs of leptons of opposite charge (opposite-sign). The search includes the ee and yu/i final 
states. 



Figure 1: The tree-level diagram for the production of a heavy Majorana neutrino {N) in the mTISM model. Lepton 
flavour is denoted by a and p. Lepton flavour is assumed to be conserved, such that a - p. The W boson produced 
from the N decay is on-shell and, in this case, decays hadronically. 

The search is guided by two theoretical models. In the first model, the SM is extended in the simplest way 
to include right-handed neutrinos [8], such that light neutrino masses are generated by a Type-I seesaw 
mechanism or by radiative corrections [9]. In this minimal Type-I seesaw mechanism (mTISM), the 
heavy Majorana neutrinos, N, can be produced via an off-shell W boson, pp (IT-)* —> i-N. Due to 
previous limits [10, 1 1], the heavy neutrino is assumed to be more massive than the IT boson and therefore 
subsequently decays to an on-shell IT boson and a lepton. The on-shell IT boson produced in the decay 
of the heavy neutrino predominantly decays into a quark-antiquark {qq) pair. The final slate in this case 
contains two opposite-sign or same-sign leptons and at least two high-px jets, where px is the transverse 
momentum with respect to the beam direction. ^ The tree-level process is illustrated in figure 1 . The free 
parameters in this model are the mixing between the heavy Majorana neutrinos and the Standard Model 
neutrinos, T^a?, and the masses of the heavy neutrinos, nif^. In this framework, LEP has set direct limits 
for wta? < mz [10, 11] and CMS has set direct limits for 90 < nij^ < 200 GeV in ee final stales [12] and 
40 < niN < 500 GeV in pp final states [13]. 

The second model is the left-right symmetric model (LRSM) [4, 14-16], where a right-handed symmetry 
SU(2 )r is added to the SM. The symmetry SU(2 )r is assumed to be the right-handed analogue of the 
SM SU(2 )l symmetry. In this model, heavy gauge bosons Tr = (ITr,Z'} are also predicted and, in this 
analysis, the heavy gauge bosons are assumed to be more massive than the heavy neutrinos, such that they 
are kinematically allowed to decay into final states including heavy neutrinos. These can be produced in 
the decays of heavy gauge bosons according to ITr —> Ni and Z' —> NN and can subsequently decay via 
an off-shell ITr boson into a lepton and a qq pair, N with IT^ ^ qq'. The tree-level processes 

are shown in figure 2. A previous ATLAS search in this framework has excluded < 2.3 TeV for 

* ATLAS uses a right-handed coordinate system, with its origin at the nominal interaction point in the centre of the detector. 
The z-axis points along the beam direction, the x-axis from the interaction point to the centre of the LHC ring, and the ^-axis 
upwards. In the transverse plane, cylindrical coordinates (r, 0) are used, where 0 is the azimuthal angle around the beam 
direction. The pseudorapidity 77 is defined via the polar angle 6 as ;/ = - In tan (0/2) . 


3 




Figure 2: The tree-level diagrams for the production of a heavy Majorana neutrino (N) in the LRSM model, in 
which heavy gauge bosons VTr and Z' are also incorporated. Lepton flavour is denoted by a and p. Lepton flavour 
is assumed to be conserved, such that a - p. The ITr boson produced from the N decay is off-shell and, in this 
case, decays hadronically. 


wiwr - mN > 0.3 TeV at 95% confidence level (CL) [17]. A more recent search performed by CMS has 
excluded nrvvR < 3.0 TeVfor > 0.05 TeV at 95% CL [18]. There are no such limits for the 

production of heavy neutrinos from Z' boson decays. 

Both the mTlSM and LRSM models produce final states containing two same-sign leptons and high-pT 
jets, but the kinematic characteristics of the events are quite different. In the mTlSM final state, one can 
reconstruct the resonant SM W boson from the jets originating from the tree-level qq pair, whereas in 
the LRSM final states, one can instead reconstruct the masses of the heavy gauge bosons. Furthermore, 
the energy scales of the two models are largely separate. The energy scale of mTlSM final states is set 
by the heavy neutrino mass, which, based on the LEP constraints [10, 11], is assumed to be greater than 
100 GeV. Instead, the energy scale of LRSM final states is set by the masses of the heavy bosons, which, 
motivated by the earlier heavy neutrino searches, are assumed to be greater than 400 GeV. For these 
reasons, the event selection criteria are optimised separately for each model, although a common object 
selection is used in both cases. 


2 The ATLAS detector 

The ATLAS detector [19] surrounds the interaction point and covers nearly the entire solid angle. The 
detector consists of an inner detector (ID) tracking system, electromagnetic and hadronic calorimeters, 
and a muon spectrometer (MS) that surrounds the other detector systems. The ID tracking system consists 
of a silicon pixel detector, a silicon microstrip tracker, both covering \q\ < 2.5, and a transition radiation 
tracker covering \q\ < 2.0. The ID tracker is immersed in a 2 T axial magnetic field provided by a 
superconducting solenoid magnet. The electromagnetic accordion calorimeter is composed of lead and 
liquid-argon (LAr) and provides coverage for \q\ < 3.2. Hadronic calorimetry is provided by steel and 
scintillator tile calorimeters for \q\ <1.7 and copper and LAr calorimeters for 1.5 < \q\ < 3.2. Additional 
LAr calorimeters with copper and tungsten absorbers cover the forward region. The MS consists of 
dedicated trigger chambers covering \q\ < 2.4 and precision tracking detectors covering \q\ < 2.7. A 
system of three superconducting toroids (one in the barrel, two in the end-caps), with eight coils each. 


4 




provides the magnetic field for the MS. A three-level trigger system is used to collect the data. The first- 
level trigger (LI) is implemented in hardware, using a subset of detector information to reduce the event 
rate to no more than 75 kHz. This is followed by two software-based trigger levels (L2 and EF), which 
together further reduce the event rate to less than 1 kHz. 


3 Background and signal simulation 

There are several SM interactions that can produce pairs of isolated charged leptons from vector bo¬ 
son decays, specifically from Drell-Yan and diboson processes and also from the decay products of top 
quarks. These processes are modelled using events produced by Monte Carlo (MC) generators. A sum¬ 
mary of the primary MC generators used in this paper to model the background processes is presented 
in table 1 . Processes that contribute to the background with pairs of same-sign leptons are indicated by 
the label ‘SS’ and processes which contribute with pairs of opposite-sign leptons, which are included in 
the charge-misidentification background estimate (section 5.2), are indicated by the label ‘OS’. The MC 
samples are normalised using the best available theoretical cross-sections, typically next-to-leading order 
or next-to-next-to-leading order in QCD. 

The production of top quark pairs and the production of a single top-quark in association with a W boson 
are simulated using MC@NLO 4.0.3 [20, 21] with the CTIO PDF set [22]. The production of single 
Z boson and diboson WZ and ZZ) events are simulated using Sherpa 1.4.1 [23] with the CTIO 

PDF set. These samples include contributions from virtual photons, with the requirement for electron 
pairs to satisfy nie+g- >0.1 GeV. No requirement is made on . The Sherpa samples include leading- 
order matrix elements for the production of up to three additional partons. The matching between the 
matrix elements and the parton shower is achieved using the CKKW method [24]. The process qq 
q'q'W-W- is generated using Madgraph 2.1.1 [25] with the MSTW 2008 PDF set [26]. The production 
of gauge bosons in association with top quark pairs is also simulated with Madgraph 2.1.1 using the 
CTEQ6F1 PDF set. 

Table 1 also shows how the various signal processes are modelled. Signal events in the mTISM model 
are generated at leading order in QCD using the Alpgen 2.14 MC generator [27, 28] with the CTEQ6F1 
PDF set. The events are generated for heavy Majorana neutrino masses between 100 and 500 GeV. Final 
states which contain exactly two prompt leptons [ ee, jiji } with same-sign charge are produced. 

The Pythia 8.170 generator [29] is used to generate FRSM events at leading order in QCD. It is assumed 
in the model that the coupling between heavy gauge bosons and the heavy neutrino is equal to the coupling 
between the respective SM gauge bosons and light neutrinos. Similarly, the couplings of the new gauge 
bosons and the quarks are assumed to be equal to the couplings between the SM gauge bosons and the 
quarks. In the Pythia implementation of this process, all of the decay products of heavy neutrinos are 
distributed isotropically and so the heavy neutrino decays are independent of my/^, with the assumption 
that my/^ > niN. In this approximation, the decay of the Z' boson is also independent of The events 
are generated with VTr boson masses between 0.6 TeV and 4.5 TeV and Z' boson masses between 0.4 TeV 
and 3.6 TeV. At each Z' and ITr mass point, the heavy neutrino mass is varied upward from 50 GeV to 
at most 100 GeV below the mass of the heavy gauge boson. At each mass point, a sample is generated 
assuming no mixing between the heavy neutrinos, which results in final states containing same-flavour 
leptons [ee, /i/r, rr}. 
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Pythia 
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CTEQ6L1 


Signal processes 


IT* ^ eN 


Alpgen 


CTEQ6L1 


SS 

Pythia 

Pythia 

MSTW 2008 

Z' NN 




Table 1; Overview of primary MC samples used for the simulation of signal and background processes. The cat¬ 
egory labelled ‘OS/SS’ refers to whether the process leads to pairs of opposite-sign (OS) or same-sign (SS) leptons. 
As described in section 5.2, OS MC samples are used in the prediction of the charge-misidentification background. 


Parton showering, fragmentation, hadronisation and the modelling of the underlying event for all Mad- 
GRAPH and Alpgen samples are performed with Pythia 8.165 and for MC@NLO samples with Herwig 
6.520 [30] and Jimmy 4.31 [31]. 

The effect of multiple pp collisions in the same or different bunch crossings is incorporated into the sim¬ 
ulation by overlaying minimum-bias events generated using Pythia 8 onto hard-scatter events, where the 
number of additional interactions is distributed in the same way as in data. All the background samples are 
produced using a simulation of the ATLAS detector [32] based on Geant4 [33]. The signal samples are 
processed through a fast simulation using a parameterisation of the performance of the ATLAS electro¬ 
magnetic and hadronic calorimeters [34], and Geant4 in the ID and MS. Both the signal and background 
samples are then processed with the same reconstruction software as the data. Small differences between 
data and MC simulation in the lepton reconstruction, identification and trigger efficiencies are corrected 
for by using specific data-driven measurements. 


4 Data sample and event selection 

The events used were selected from pp collision data with an integrated luminosity of 20.3 fb“^ collected 
by ATLAS in 2012. Quality criteria are applied to suppress non-collision backgrounds such as cosmic-ray 
muons, beam-related backgrounds, and spurious noise in the calorimeters. 


4.1 Object reconstruction and selection 

The search uses reconstructed electrons, muons, jets and a measurement of the missing transverse mo¬ 
mentum. 
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Electrons are required to satisfy tight identification requirements [35] and to have pj > 20 GeV and 
It/I < 2.47. Any electron in the transition region between the barrel and end-cap calorimeters (1.37 < \r]\ < 
1.52) is rejected. In order to avoid double counting electrons as jets, the nearest jet within A/?(e, jet) = 0.2 
of an electron and with pj < lEj, where Ej - E sin 6 is the transverse energy deposited by the electron, 
is rejected. 

Muons are required to be reconstructed in the MS and successfully matched to a good-quality track in 
the ID [36]. It is required that muons have pj > 20 GeV and \r]\ < 2.5. In order to suppress muons with 
misidentified charge, it is required that there is a consistent measurement of charge in the MS and ID. 
Muons with pj < 80 GeV are required to be well separated from jets, such that A/?(yu,jet) > 0.4, where 
AR - + (A0)2. 

Jets are reconstructed using the anti-kf clustering algorithm [37, 38] with the radius parameter set to 0.4. 
Jets are calibrated [39, 40] using an energy- and Ty-dependent simulation-based calibration scheme, with 
in-situ corrections based on data. The impact of multiple overlapping pp interactions is accounted for 
using a technique that provides an event-by-event and jet-by-jet correction [41]. Events are rejected if 
any jet is identified as originating from beam-halo effects or calorimeter noise. Jets are required to have 
Pi > 20 GeV and \t]\ < 2.8. The pi requirement is chosen in order to maximise the acceptance for the 
mTISM model. Eor jets with pi < 50 GeV within the acceptance of the tracking detector (\t]\ < 2.4), the 
‘jet vertex fraction’ (JVE) [42] is required to be greater than 0.5. The jet vertex fraction is calculated by 
summing the pi of tracks associated with the jet and matched to the selected primary vertex, and dividing 
it by the sum of the pi of all tracks associated with the jet. 

The primary vertex of the event is defined as the reconstructed vertex with the highest X P^, consistent 
with the beam spot position, where the sum is over all tracks associated with the candidate primary 
vertex. 

The missing transverse momentum, is used to identify invisible particles such as light neutrinos 

that escape detection. The E™®® quantity is calculated as the magnitude of the negative vector sum of 
all reconstructed particles momenta, including muons, electrons, photons, and jets, as well as clusters of 
calorimeter cells, not associated with these objects. 


4.2 Lepton isolation criteria 

Backgrounds due to misidentified leptons and non-prompt leptons, which are described in detail in sec¬ 
tion 5.3, can be suppressed by requiring that leptons are isolated from other activity in the event. Because 
of the different background compositions for electrons compared to muons and of the different response 
of the detector to isolated electrons and muons from prompt sources, different isolation criteria are used 
for the two lepton flavours. 

Electrons are required to satisfy (^p^^ + 1 GeV^ /Ej < 0.05 and E^^ IEj < 0.05, where p^^ is the sum of 
the Pi of all tracks within a cone of AR = 0.2 around the electron, excluding the electron track itself, and 
Ej^ is the sum of the transverse energy in a cone of AR = 0.3 around the electron, excluding the electron 
itself. The criteria are looser at high electron Ei in order to maintain high efficiency. 

Muons with pi < 80 GeV are required to have p^^Ipi < 0.05 and E^^lpi < 0.05, where p^^ is the 
sum of the pi of all tracks within a cone of AR = 0.3 around the muon, excluding the muon track itself, 
and E^^ is the sum of transverse energy measured in the calorimeter within a cone of AR = 0.2 around 
the muon, excluding energy deposits associated with the muon track. Eor muons with pi > 80 GeV 
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the requirements are relaxed in order to maintain high efficiency. Muons with pj > 80 GeV are either 
required to satisfy pj < 0.05, or if they are within AR = 0.4 of a jet they can additionally be selected 
if {mfjj - rtij) > 10 GeV, where rtij is the reconstructed mass of the jet closest to the muon and m^j is the 
invariant mass of the jet and the muon. The latter criterion is efficient for the decay of a boosted heavy 
neutrino decaying into a muon and a qq pair, while rejecting a large fraction of misidentified muons. 

Both the muons and electrons must satisfy a set of requirements on the impact parameters at the primary 
vertex in order to further suppress leptons originating from heavy-flavour decays. They are required to 
have a transverse impact parameter, d^, which satisfies |r/ol < 0.2 mm and |r/ol/o‘(‘5?o) < 3, where cr(r/o) is 
the uncertainty on d^. It is also required that the product of the longitudinal impact parameter (zq) and the 
sine of the polar angle of the lepton {9) satisfy |zo sin 6\ <2 mm. 


4.3 General event selection 

The events are required to satisfy one of a suite of triggers [43] that select events with either one or two 
high-pT leptons. This analysis uses single-lepton {e or p) triggers with a 24 GeV px threshold. The 
analysis also uses a dimuon trigger for events in which one muon has satisfied a px threshold of 20 GeV 
and a second muon has satisfied a threshold of 8 GeV. The choice of triggers is found to maintain the 
highest possible signal efficiency in each channel across the presented range of heavy neutrino masses. 
The electron (muon) trigger efficiencies for offline selected electrons (muons) are > 94% (70%) and > 
85% (90%) in the barrel and end-cap, respectively. The total efficiency for a single electron (muon) within 
the detector acceptance to satisfy the full lepton selection described in sections 4.1 and 4.2, including the 
trigger requirement, is approximately 54% (70%). 

The highest-px lepton in an event is required to satisfy px > 25 GeV. The choice of lepton px threshold 
is dictated by the trigger requirements. It is required that at least one lepton with px > 25 GeV is matched 
to one of the described triggers. In the case of the event being selected by the dimuon trigger, two muons 
must be matched to the trigger. Any other leptons must satisfy px > 20 GeV. 

Events are required to contain exactly two leptons from {ee, pp) with same-sign charge, where the two 
leptons must have ID tracks associated with the same vertex. To remove the small background arising 
from muon bremsstrahlung in the ID or in the first layers of the EM calorimeter, events are rejected 
if a muon’s ID track is also reconstructed as an electron. Backgrounds from ITZ and ZZ decays are 
suppressed by rejecting events which contain an additional lepton, where the additional lepton is selected 
with looser identification requirements and no requirements on the isolation variables. The impact of 
the latter criterion on the signal efficiency is negligible compared to the overall uncertainty on the signal 
acceptance. 


4.4 Selection criteria for mTISM signal events 

The signal region for the mTISM model is defined for events containing, in addition to the two leptons, 
at least two jets. The invariant mass, mp, of the two highest-px jets is required to lie in the range 
60 < nijj < 100 GeV. This selects events consistent with an on-shell W boson decaying to a qq pair. 
The invariant mass of the two leptons {mu) is required to be greater than 40 GeV, a selection that has 
high efficiency for the signal over the full range of values considered. In the ee channel, the charge- 
misidentification background (described in section 5.2) is suppressed by requiring that the invariant mass 
of the two leptons is outside a window around the Z boson mass, \mu - mz\ > 20 GeV. Backgrounds 


due to electroweak processes producing same-sign leptons are dominated by those including at least one 
light neutrino, particularly those arising from diboson production processes. Such backgrounds may 
have high compared to the mTISM signal, so the events in the signal region are required to have 
< 40 GeV. 

The total efficiency (including the detector acceptance) for signal events to satisfy all selection criteria 
is lower in ee events than in ///r events as electrons have lower efficiency fo safisfy the identification 
and isolation criteria. The efficiency for leptons in mTISM signal events to satisfy the object selection 
increases as a function of the lepton pj. The total efficiency therefore increases as a function of from 
approximately 0.5% to 24% in the ee channel and 3% to 30% in the pp channel. 


4.5 Selection criteria for LRSM Wr and Z' signal events 

The signal region for heavy neutrinos produced in the decays of ITr bosons is defined for events contain¬ 
ing at least one jet. To exploit the high energy scale of the signal events it is required that mu >110 GeV 
and the invariant mass of the system consisting of the two leptons and one or two jets must satisfy 
^UjU) > 400 GeV. If the event contains more than two jets, the two highest-px jets are used. 

The signal region for heavy neutrinos produced in the decays of Z' bosons uses the same requirement on 
the dilepton mass, but requires at least two jets and the invariant mass of the system consisting of the two 
leptons and two to four jets must satisfy mujj(jj) > 200 GeV. If the event contains more than four jets, 
the four highest-px jets are used. A lower invariant mass of the system of 200 GeV is considered in the 
search for Z' bosons compared to the value of 400 GeV used for ITr bosons since previous searches for 
heavy neutrinos have already set strong constraints for otwr < 400 GeV [17], whereas there are no such 
limits for Z' production. 

Allowing only one jet in the VTr decay or two jets in the Z' decay (rather than two and four respectively) 
increases the signal efficiency for the case when the heavy neutrino is boosted (mvg^ » ni/v) so that the 
qq pair produced in the decay of the off-shell Wr results in a single jet reconstructed in the detector. For 
niy^ » mu, up to 60% of events contain only one jet in ITr events or two jets in the case of Z' decays. The 
total selection efficiency for LRSM events depends on the masses my^ and m^ and also the ratio my^lmu- 
The total efficiency (including the detector acceptance) ranges from approximately 0.5% to 25% in the 
ee channel, and from approximately 1.5% to 30% in the pp channel. For small values of m^lmy^, the 
efficiency decreases rapidly and is below 15% for mulmvg^ < 0.1, because the heavy neutrino decay 
products are highly boosted and the leptons are less isolated. 


5 Background estimation 


The background is evaluated in the signal regions according to three categories of lepton pairs. The 
first two categories describe the contribution due to leptons originating from SM processes that produce 
prompt isolated leptons. The first category, labelled as ‘prompt’ in the following and discussed in sec¬ 
tion 5.1, corresponds to irreducible background from true same-sign prompt lepton pairs. The second cat¬ 
egory, labelled as ‘charge-flip’ in the following and discussed in section 5.2, corresponds to true opposite- 
sign prompt lepton pairs, in which one lepton has its charge mismeasured. The third category, labelled 
as ‘non-prompt’ in the following and discussed in section 5.3, corresponds to one or both leptons being 
either a non-prompt lepton from semileptonic heavy-flavour decays, a jet misidentified as a lepton or, in 
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the case of electrons, photons misidentified as leptons. The first category of background events are en¬ 
tirely estimated using the SS MC samples described in section 3, the second category is estimated using 
OS MC samples in conjunction with a measurement in the data of the charge misidentification rate and 
the third category is estimated from data. 


5.1 Background from prompt same-sign leptons 

The background from SM processes that lead to two same-sign prompt leptons is referred to as the prompt 
background and is estimated using the MC samples described in section 3. The largest contribution to 
the SM background in the signal regions originates from ITZ and ZZ events. Other prompt contributions 
were estimated to be negligible; these include processes involving the production of three electroweak 
gauge bosons or of a Higgs boson. The simulation of the ITZ and ZZ backgrounds is validated by select¬ 
ing events with either three or four charged leptons satisfying the selection cuts described in section 4.3. 
These events arise predominantly from ITZ and ZZ production, respectively, with a negligible contribu¬ 
tion expected from the signal processes. The expected and observed jet multiplicity and leading jet pj 
distributions for these events are shown in figure 3. The number of events predicted by the simulation is 
found to be in good agreement with the data and the kinematic properties of the events are adequately 
described by the simulation. 



1 2 3 4 >5 

Number of jets 



(a) 


(b) 


Figure 3: Distribution of (a) the number of jets and (b) the leading jet pj in events containing any combination 
of exactly three or four leptons. The events must contain one lepton with > 25 GeV and all other leptons 
must satisfy pj > 20 GeV. The contribution labelled ‘Other’ is from processes described in section 5.1 (with MC 
samples described in section 3), other than the contributions from WZ and ZZ, which are labelled separately. The 
shaded bands indicate the experimental uncertainties on the total expected background, including all contributions 
described in section 6.2, but not including any uncertainty on the WZ and ZZ cross-sections. The lower plots show 
the ratio of data to the total expected background. 
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5.2 Background from prompt opposite-sign leptons 

SM processes that produce opposite-sign leptons can also enter into the signal regions if the charge of 
one lepton is incorrectly measured in the detector. This is referred to as the ‘charge-flip’ background and 
includes pairs of opposite sign-leptons produced in tt, W-W'^ and Z processes. The MC samples used to 
model these processes are described in section 3. 

The probability for a charge-flip event to occur is measured using Z ^ fifi and Z —> cc events, which are 
identified using events, which have exactly two leptons with an invariant mass close to the mass of the Z 
boson. For muons, the charge-flip rate is estimated by comparing the two independent measurements of 
the muon charge in the MS and the ID. The charge-flip rate is found, as expected from simulation, to be 
consistent with zero, and contributes a negligible number of events to the expected background. 

There is however a sizeable charge-flip rate for electrons, due to bremsstrahlung photons produced in the 
ID and subsequently converting to electron pairs. A sample of Z —> ce events in data is used to perform 
a single maximum-likelihood fit, which extracts the electron charge-flip rate as a function of rj [44]. The 
contribution to these events from fake and non-prompt electrons (see section 5.3) is subtracted from the 
data sample prior to the fit. The measured rate, which is strongly correlated to the amount of material 
traversed by the electron in the ID, is found to be approximately I0“‘^ for electrons within the barrel region 
(l'7l < 1.0), increasing to 10 ^ at the edge of the detector acceptance (I 77 I = 2.47). The data-to-MC ratio 
of the measured electron charge-flip rate is used as a correction factor to the charge-flip rate obtained 
in the simulation. The pj-dependence of the charge-flip rate is therefore directly taken from the MC 
simulation. 


5.3 Background from fake and non-prompt leptons 

Events where jets or photons are misidentified as leptons (‘fakes’) or events with non-prompt leptons 
which originate from semileptonic heavy-flavour decays constitute a significant background, which is 
referred to as the non-prompt background. These processes include W + jets and tt production, where 
one lepton originates from a vector boson decay and the other lepton is misidentified or from a non-prompt 
decay. This background cannot be reliably predicted from MC simulation and is estimated directly in each 
of the signal regions, using the data-driven matrix method [45]. The matrix-method characterises leptons 
which satisfy ‘loose’ identification criteria as being from prompt or fake/non-prompt sources according 
to their probabilities to satisfy the full lepton identification criteria. 
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A ‘loose’ electron is defined with an identification requirement that is relaxed from the ‘tight’ to ‘medium’ 
operating point compared to the standard electron selection described in section 4.1. The selection criteria 
for a ‘loose’ muon are identical to the full selection described in section 4.1, with the same requirements 
on the impact parameters described in section 4.2 but without a requirement to satisfy any isolation 
criteria. 

The measurement of the probabilities r and / for ‘loose’ leptons from prompt or fake/non-prompt sources 
respectively to satisfy the full lepton identification criteria are the key factors in the estimate of the non¬ 
prompt backgrounds. The probabilities r are measured in Z ^ events as these events are dominated 
by prompt leptons. 

The probabilities / are measured in a selection of events which contain a large number of fake and 
non-prompt leptons. In these control samples, any residual prompt background is subtracted using MC 
estimates, and low-mass hadronic resonances are excluded by requiring niff >15 GeV. For electrons, 
the probabilities are measured in events with at least one jet and with exactly one electron. To suppress 
events containing W decays it is required that |A0(e, < 0.5 and + mT < 40 GeV, where nij 

is the transverse mass.^ The probabilities for muons are measured in events containing pairs of muons 
chosen with transverse impact parameter requirements which differ from the standard selection, such that 
|r/ol < 10 nim and \do/cr{do)\ > 5. This sample of pairs of muons is expected to have a composition similar 
to the signal region, making it suitable for the measurement of /. A correction factor is applied to / for 
muons to account for the fact that muons with high \do/(T{do)\ have a lower probability than muons with 
low \do/cr{do)\ to satisfy the isolation criteria. The correction factor of approximately 1.4 is measured 
using bb and tt MC simulation. 

A validation region is defined for same-sign nn events which have satisfied the dimuon trigger (section 
4.3) with one muon satisfying pj > 25 GeV and the second satisfying 10 GeV < pT < 20 GeV. As the pr 
of the second muon is relatively low, this region is dominated by non-prompt background events. The jet 
multiplicity measured in events in this region is compared to the expected background as shown in figure 
4. The overall level of agreement is within one standard deviation for up to five jets in the event. 

5.4 Validation of background estimates 

The validation of prompt, charge-flip and non-prompt background estimates is considered in this sec¬ 
tion. 

The combined background estimate can be evaluated using events containing exactly two same-sign 
leptons and no jets. This sample of events is orthogonal to each of the different signal regions described in 
sections 4.4 and 4.5, and is expected to contain only a negligible contribution from possible signal events. 
Comparisons of the distributions of the F"™®* and lepton pj as measured in data and estimated from the 
background predictions from the three sources described above are shown in figure 5. The overall agree¬ 
ment between data and background prediction is within approximately one standard deviation in both the 
ee and pp channels. 

An additional sample of events is considered, which fully includes all signal regions. This region is 
defined for events containing exactly two same-sign leptons, with no requirement on the number of jets in 
the event. The contamination from signal events in this sample is less than 2% of the total background in 


■ iri'Y = 


^2pT£:r'(i - cos [e, where pj refers to the electron. 
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Figure 4; Distribution of the number of jets in a validation region consisting of events containing exactly two same- 
sign muons with one muon satisfying pi > 25 GeV and the second satisfying 10 GeV < px < 20 GeV. The 
shaded bands indicate the total uncertainty, including all contributions described in section 6, on the total expected 
background and the lower plots show the ratio of data to the total expected background. 


both channels. The distributions of the leading jet pj and the number of jets in these events, in the ee and 
pp channels, as measured in data and estimated from the background predictions from the three sources 
described above, are shown are shown in figure 6. The overall agreement between data and prediction in 
the two channels is within approximately one standard deviations. 
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Figure 5; The (top) and lepton pj (bottom) distributions for the ee (left) and pp (right) channels in a validation 
region consisting of events with exactly two same-sign leptons and no jets. The shaded bands indicate the total 
uncertainty, including all contributions described in section 6, on the total expected background and the lower plots 
show the ratio of data to the total expected background. 
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Figure 6; The distribution of the transverse momentum pj of the leading jet (top) and the distribution of the number 
of jets (bottom) for the ee (left) and pp (right) channels in a validation region consisting of events with exactly two 
same-sign leptons. The shaded bands indicate the total uncertainty, including all contributions described in section 
6, on the total expected background and the lower plots show the ratio of data to the total expected background. 
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6 Systematic uncertainties 


The background estimates and signal efficiencies are subject to several systematic uncertainties. The 
relative size of the uncertainties on the total background estimates in the mTISM and LRSM signal regions 
are detailed in tables 2 and 3 respectively. Since the tables show the size of each uncertainty relative to 
the total background estimate, the impact of each uncertainty depends on the background composition in 
the different channels (see section 7 and tables 4 and 5). 


6.1 Background uncertainties 

The systematic uncertainty on the estimate of the non-prompt background is dominated by the uncertain¬ 
ties on the measurements of the rates for leptons from fake and non-prompt sources to satisfy the lepton 
identification criteria. This uncertainty is dominated by the effect of the choice of control region definition 
on the measured probabilities, as well as the statistical uncertainty on determining the probabilities. The 
total uncertainty on the non-prompt background in the ee (jj.fi) channel varies from 30 - 48% (41 - 45%), 
depending on the signal region. 

The charge-flip background is only relevant for the ee channel (see section 5.2) and its uncertainty is 
dominated by the statistical precision with which the charge-flip rate is determined from the available 
data. There is, additionally, a non-negligible contribution to the total uncertainty due to the modelling of 
the subtraction of the non-prompt background, although this is correlated with the uncertainty on the non¬ 
prompt background. Since the charge-flip background uses MC simulafion, the systematic uncertainties 
discussed in section 6.2 also apply to the charge-flip background. The fotal uncertainty on the charge-flip 
background varies from 18 - 46%, depending on the signal region. 

The uncertainty on the normalisation of the backgrounds originating from WZ and ZZ processes is de¬ 
rived from the diboson control region described in section 5.1. This uncertainty is taken to be either the 
difference between the data and the prediction or the statistical uncertainty from the limited data stat¬ 
istics, whichever is largest. The uncertainty is applied as a function of the number of jets, leading to 
uncertainties of 10 - 14%, depending on the signal region. The uncertainty on the cross-section of other 
background processes from MC estimates is taken from the uncertainties on the theoretical cross-sections. 
The combined effect of the diboson normalisation and theoretical cross-section uncertainties is labelled 
as ‘Prompt normalisation’ in tables 2 and 3. 

Another source of uncertainty is due to the MC statistical uncertainty. In the mTISM signal region this is 
particularly large in the ee channel due to the small number of events in the Z ^ ee MC sample in this 
region. 


6.2 Uncertainties on MC simulation 

The following uncertainties are applied to all MC-derived predictions. In addition to affecting the prompt 
and charge-flip background esfimates, these uncertainties also apply to the signal simulation. 

The systematic uncertainty on the jet energy scale has an important effect on both the signal and back¬ 
ground processes, as the signal-region event selection includes requirements on quantities that are recon¬ 
structed from jet kinematics [40]. The uncertainty on the jet energy scale is notably asymmetric in the 
mTISM signal region (table 2) compared to the LRSM signal region (table 3) due to the dependency on 
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ee 

nn 

Non-prompt 

±1 

+ 14 

Charge-flip 

±1 

- 

Prompt normalisation 

± 2 

+ 10 

MC statistics 

+ 14 

+ 8 

Jet energy scale 

H-8/-22 

-H7/-12 

^miss 

+21-3 

H-6/-3 

Jet energy resolution 

+ 4 

+ 3 

Jet vertex fraction 

+21-6 

H-4/-5 

Lepton uncertainties 

± 2 

+21-3 

Luminosity 

± 2 

± 1 

Total 

+201-29 

+211-24 


Table 2: A breakdown of the relative uncertainty on the total background (given in %) in the mTISM signal region. 
The various sources of systematic uncertainty are described in Section 6. 



Wr 

ee 

Z' 

ITr 

Z' 

Non-prompt 

+ 12 

± 15 

+ 10/-9 

+ 13/-11 

Charge-flip 

-H5/-4 

+41-3 

- 

- 

Prompt normalisation 

+ 5 

+51-4 

± 18 

+ 14 

MC statistics 

+ 7 

+ 6 

±5 

+ 4 

Jet energy scale 

+81-9 

+51-4 

+ 5 

+5/-4 

Jet energy resolution 

±0.9 

+ 0.6 

± 1.2 

+ 0.2 

Jet vertex fraction 

+ 1 

+ 2 

+0.9/-0.1 

+2.4/-1.1 

Lepton uncertainties 

+21-1 

+2.7/-1 

± 2 

± 2 

Luminosity 

± 1 

+ L5/-1 

± 1 

+1 

Total 

+ 18 

± 18 

+ 21 

+20/-19 


Table 3; A breakdown of the relative uncertainty on the total background (given in %) in the LRSM signal regions. 
The various sources of systematic uncertainty are described in Section 6. 
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the leading dijet mass in the mTISM signal region definition described in section 4.4. The uncertainty 
due to the JES on the sum of the prompt and charge-flip backgrounds varies from 6 - 28%, depending on 
the signal region. The modelling of the missing transverse momentum (‘E™®®’) [46] is included as a sys¬ 
tematic uncertainty in the mTISM signal region (table 2) as there is a dependence on this quantity in the 
signal region definition described in section 4.4. Other smaller systematic uncertainties include the un¬ 
certainty on the the jet energy resolution [47], the jet vertex fraction requirement, uncertainties on lepton 
identification efficiencies, energy / momentum scales and resolutions (‘Lepton uncertainties’) [35, 36] 
and the uncertainty on the luminosity measurement (‘Luminosity’) [48]. 


6.3 Signal-specific modelling uncertainties 

In addition to the uncertainties associated with the MC simulation of background processes, there are 
systematic modelling uncertainties associated with the signal MC samples. An uncertainty is considered 
for the signal MC simulation to reflect the choice of parton shower model. The nominal parton shower 
model that is used for all signal MC samples is Pythia 8.165. The total number of events in the signal 
region when the signal MC generator is interfaced to Pythia is compared to the number of events when 
the generator is interfaced to Herwig 6.520. The variation in the signal efficiency is measured to be 
approximately 5%. The uncertainty due to the parton distribution functions on the signal acceptance is 
found to be approximately 5% for the mTISM signal samples and approximately 7% for the LRSM signal 
samples. 

A systematic uncertainty is also considered to cover the effect of using the fast detector simulation de¬ 
scribed in section 3. Two versions of the MC signal, one with a full detector simulation and the other with 
the standard fast detector simulation, are compared in each lepton channel, for a single signal mass point. 
The difference in the efficiency to select signal events, approximately 4%, is assigned as an uncertainty. 


7 Results 

The numbers of events measured in data are compared to the expected numbers of background events 
in the signal regions, with the intention of interpreting an excess of events in data in terms of a heavy 
Majorana neutrino in the mTISM or LRSM models. 


7.1 Results in the mTISM signal region 

The observed and predicted distributions of the invariant mass of the two highest-pT jets {mjj) in events 
satisfying the mTISM selection criteria, excluding the criteria on rujj (60 GeV< rujj < 100 GeV) are 
presented in figure 7. The shapes of the distributions show good agreement between data and expectation. 
The numbers of expected and observed events in the mTISM signal region (indicated by arrows in figure 
7) are shown in table 4. There is no excess of events relative to the expectation. The observed yields 
in the data are used to set 95% CL upper limits on the production cross-section times branching ratio, 
cr X Br(pp — > i-N — > £-i-qq'), of heavy neutrinos to electrons or muons, using the profile-likelihood 
test statistic [49, 50] and the CLs method [51]. The systematic uncertainties are included in the test 
statistic as nuisance parameters. Each systematic uncertainty is assumed to be uncorrelated with all other 
systematic uncertainties. The limits are shown as a function of mi^ in figure 8 and are translated into limits 
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on the mixing parameter, \VtN?', between the heavy neutrino and the SM neutrino, separately for the ee 
and njj. channels. The extraction of the limits on the mixing parameters uses the leading-order cross- 
section for the signal process and no uncertainties are included on the signal cross-section. The extraction 
furthermore assumes that only the lightest of the heavy neutrinos contributes to the cross-section and 
that the masses of the other heavy neutrino species are sufficiently high that the effect of interference is 
negligible. 
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Figure 7: Invariant mass of the two highest-px jets (nijj) in events satisfying the mTISM signal region criteria 
(excluding the mjj criteria) for (a) ee and (b) mi events. Events satisfying all selection criteria are in the region 
indicated by the arrows. The expected mTISM signal disttibution for - 120 GeV is shown by the dashed 
(blue) histogram. The values of the mixing parameter are chosen such that the signal distribution is clearly 
visible, \VeN?' - 0.2 for (a) and IVjuArp = 0.02 for (b). The shaded bands indicate the total uncertainty, including all 
contributions described in section 6, on the total expected background and the lower plots show the ratio of data to 
the total expected background. 



ee 

pp 

Prompt 

o c-fO.9 

5.8!1;3 

Charge-flip 


< 0.02 

Non-prompt 

4.3 + 1.8 

2.9+ 1.3 

Total background 

21+"^ 

8.7+ 2.0 

Data 

19 

6 

Signal (m^v = 120 GeV) 

18 

6.6 

Signal {niN - 240 GeV) 

30 

5.3 


Table 4: Total event yields measured in data and predicted for signal and background processes in the mTISM signal 
region. The uncertainties shown on the various backgrounds correspond to the total uncertainty. The expected 
number of mTISM signal events are calculated for = 120 GeV (with IVmP equal to 0.03 and 0.003 in the ee 
and pp channels respectively) and nit^ - 240 GeV (with |V7ivP equal to 0.2 and 0.02 in the ee and pp channels 
respectively). The values of the mixing parameters \VtN?' are chosen to be close to the expected limit shown in 
figure 8. For backgrounds that have zero expected events, the 95% CL upper limit is shown. 
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(c) (d) 

Figure 8: Observed and expected 95% confidence level upper limits on the cross-section times branching ratio for 
the production of mTISM heavy Majorana neutrinos as a function of the heavy neutrino mass for (a) the ee channel 
and (c) the pp channel. The limits on the mixing between the heavy Majorana neutrinos and the SM neutrinos are 
shown in (b) and (d). Values larger than the solid black line are excluded by this analysis. 


7.2 Results in the LRSM signal region 

The observed and expected numbers of events for the LRSM signal regions are shown in table 5. There 
are no excesses observed above the expected numbers of background events. 

The LRSM signal is expected to produce a peak in the invariant mass of the decay products of the heavy 
gauge boson. This would be observed in the invariant mass distribution in the ILr (Z') 

signal regions, as described in section 4. The observed and predicted distributions are shown in figures 9 
and 10. Binned likelihood fits are performed to the invariant mass distributions and the profile-likelihood 
fesf sfafisfic is used fo assess fhe compafibilify of fhe dafa wifh fhe background-only and signal-plus- 
background hypofheses. No significanl excess is observed in fhe dafa compared fo fhe background ex- 
pecfafion and 95% CL upper limifs on fhe cross-secfion of fhe producfion of heavy gauge bosons decaying 
fo heavy neufrinos wifhin fhe LRSM are sef using fhe CLg mefhod. The expecfed and observed cross- 
secfion exclusion limifs as a funcfion of fhe masses of fhe heavy gauge bosons and heavy neufrino are 
shown for example mass poinfs for bofh channels, ee and /r/r, in fable 6. The full cross-secfion limifs 
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for all analysed mass points are available in HepData.^ Exclusion contours are also set in the {my^, mi^] 
plane, and are shown in figure 11 . The contours are found by comparing the cross-section limits to the 
leading-order cross-sections for the signal processes and no uncertainties are included on the signal cross- 
sections. For this interpretation there is assumed to be no mixing between lepton flavours and three cases 
are investigated: the first two cases assume a single heavy neutrino being kinematically accessible and 
being of either electron or muon flavour, which leads to events expected in one of the ee and fijj. channels, 
respectively. The limits are stronger in the /r// channel at low heavy neutrino mass, due to the higher 
signal acceptance discussed in Section 4.5. As the efficiency of the LRSM signal is low in the ee channel 
for my^ » niN, the limit in this region is weaker. Finally the case where two degenerate neutrinos are 
present is investigated, leading to events expected in both the ee and //// channels. The expected limits for 
this scenario are slightly stronger than either of the individual channel scenarios. The sharp changes in 
the expected and observed limits visible in figure 11 originate from the interpolation between the limited 
number of mass points for which the limits are extracted. 
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Wr 

Z' 

ITr 

Z' 

Prompt 

26 + 5 

34 + 6 

33 + 8 

42+ 10 

Charge-flip 

44+11 

44+10 

^-8 

< 0.03 

<0.03 

Non-prompt 

23 + 11 

OO-Pll 

9.8!^ 

ITJ, 

Total background 

93 + 16 

111+16 

43 + 9 

uu i2 

Data 

94 

106 

44 

55 

Signal 

2.4 

5.2 

2.9 

7.7 


Table 5: Total event yields measured in data and predicted for signal and background processes in the LRSM signal 
regions. The uncertainties shown on the various backgrounds correspond to the total uncertainty. The number of 
LRSM signal events are calculated for ota') = {2600; 1950} GeV and {otz-; m^} — (2200; 550} GeV. For 
backgrounds that have zero expected events, the 95% CL upper limit is shown. 


^ http://hepdata.cedar.ac.uk/ 
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(a) (b) 

Figure 9; Invariant mass of two leading leptons and up to two leading jets after applying additional VFr selection 
criteria (two same-sign leptons, at least one jet, nia >110 GeV and muj(j) > 400 GeV), for the ee-channel (a) and 
p/t-channel (b). A finely binned LRSM signal sample is represented by the dashed (blue) histogram corresponding 
to = 2600 GeV and = 1950 GeV. The shaded bands indicate the total uncertainty, including all contribu¬ 
tions described in section 6, on the total expected background and the lower plots show the ratio of data to the total 
expected background. 




[GeV] 


(a) (b) 

Figure 10: Invariant mass of two leading leptons and up to four leading jets after applying additional Z' selection 
criteria (two same-sign leptons, at least two jets, nia >110 GeV and mujjijj) > 200 GeV), for the ee-channel (a) and 
pp-channel (b). A finely binned LRSM signal sample is represented by the dashed (blue) histogram corresponding 
to niz’ - 2200 GeV and = 550 GeV. The shaded bands indicate the total uncertainty, including all contributions 
described in section 6, on the total expected background and the lower plots show the ratio of data to the total 
expected background. 
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mv^ 

[GeV] 

niN 

[GeV] 

Efficiency [%] 

Expected 
limit [fb] 

Expected 
limit ±lcr [fb] 

Expected 
limit ±2(t [fb] 

Observed 

limit [fb] 

VTr ee 

1400 

50 

0.80 

100 

130/70 

190/49 

65 

1400 

1050 

22 

2.2 

3.1/1.5 

3.8/1.1 

1.7 

2300 

50 

0.39 

400 

530/250 

780/160 

320 

2300 

1725 

23 

0.87 

1.1/0.64 

1.8/0.44 

0.56 

2900 

50 

1.2 

140 

210/98 

270/71 

160 

2900 

2175 

23 

0.85 

1.1/0.62 

1.7/0.43 

0.54 


1400 

50 

19 

1.8 

2 .8/1.2 

3.5/0.88 

2.1 

1400 

1050 

26 

1.2 

1.7/0.83 

2.2/0.59 

1.1 

2300 

50 

18 

1.4 

2.2/0.91 

3.2/0.63 

2.1 

2300 

1725 

26 

0.85 

1.3/0.55 

1.8/0.38 

1.3 

2900 

50 

15 

2.2 

3.5/1.4 

5.4/0.99 

3.4 

2900 

2175 

25 

0.87 

1.4/0.57 

1.8/0.39 

1.3 

Z' ee 

1000 

50 

0.27 

500 

740/340 

920/240 

440 

1000 

375 

14 

5.2 

6.6/3.8 

8.3/2.8 

4.2 

2000 

50 

0.091 

1100 

1600/770 

2000/560 

1300 

2000 

750 

17 

2.8 

3.6/2.0 

4.8/1.5 

2.6 

2400 

50 

0.16 

570 

780/400 

970/290 

550 

2400 

900 

17 

2.9 

3.6/2.0 

4.8/1.5 

2.6 


1000 

50 

11 

4.8 

7.2/3.2 

10 /2.2 

8.5 

1000 

375 

24 

2.2 

3.2/1.5 

4.6/1.1 

3.8 

2000 

50 

15 

2.7 

4.0/1.8 

5.8/1.3 

3.4 

2000 

750 

23 

1.7 

2 .6/1.2 

3.3/0.83 

2.2 

2400 

50 

12 

3.7 

5.6/2.5 

7.0/1.8 

4.8 

2400 

900 

24 

1.7 

2.5/1.2 

3.3/0.82 

2.2 


Table 6: Expected and observed 95% CL limits for the production of LRSM VTr or Z', with subsequent decay into 
a heavy Majorana neutrino in the ee and pyu channels for several example mass points. Also shown is the event 
selection efficiency (including the detector acceptance). The full results for all analysed mass points are available 
in HepData. 
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Figure 11: Observed and expected exclusion contour at 95% confidence level as a function of the mass of a heavy 
Majorana neutrino and of a IFr (left) or Z' boson (right) within the LRSM. The limits in (a) and (b) show the 
scenario where the heavy neutrino has electron flavour and those in (c) and (d) show the scenario where it has muon 
flavour. The limits in (e) and (f) show the case of two degenerate neutrinos, one has electron flavour, and the other 
muon flavour (no mixing between lepton flavours is assumed). 
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8 Conclusions 


The proton-proton collision data sample with a centre-of-mass energy of 8 TeV collected by ATLAS 
with an integrated luminosity of 20.3 fb“^ has been used to search for the production of heavy Major- 
ana neutrinos. The selected events contain two same-sign charged leptons and high-px jets. Two final 
selection criteria were used to provide sensitivity to two different benchmark models, the first being a 
minimal extension of the SM, mTISM, and the second being a left-right symmetric extension of the SM, 
LRSM. The background expectation in the signal regions is estimated using a combination of data-driven 
methods and MC simulation. 

No significant excess of events is observed above the SM background, and limits are set for each model. 
In the framework of a mTISM, limits are set on the cross-section times branching ratio and the mixing 
between heavy neutrinos and SM neutrinos as a function of the heavy neutrino mass. The observed limits 
exclude mixing parameters as low as = 0.0028 and \VeN?' - 0.029. In the framework of an LRSM, 
limits are set on heavy Majorana neutrino masses in the range 50 GeV to 2000 GeV for heavy gauge boson 
masses above 400 GeV. 
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